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Abstract 

We show that the five- dimensional Maxwell theory with the Chern-Simons term is tachyonic 
in the presence of a constant electric field. When coupled to gravity, a sufficiently large 
Chern-Simons coupling causes instability of the Reissner-Nordstrom black holes in anti-de 
Sitter space. The instability happens only at non-vanishing momenta, suggesting a spatially 
modulated phase in the holographically dual quantum field theory in (3 + 1) dimensions, with 
spontaneous current generation in a helical configuration. The three-charge extremal black 
hole in the type IIB superstring theory on AdS 5 x S* 5 barely satisfies the stability condition. 



1 Introduction 



Instability of black holes in anti-de Sitter space has attracted much attention recently due to 
its relevance to quantum phase transitions in dual strongly interacting quantum field theory 
at finite density [TJ [2] . In this paper, we point out a new type of instability caused by the 
Chern-Simons term. A novel feature is that the instability happens only at non-vanishing 
momenta, suggesting a spatially modulated phase transition in the holographically dual field 
theory. 

In three dimensions, the Maxwell theory becomes massive when the Chern-Simons term 
is included [HI H] . In higher dimensions, the Chern-Simons term starts with a higher power in 
gauge fields, but it can contribute to quadratic fluctuations if there is a non-zero background 
gauge field. In this paper, we will show that the Maxwell theory in five dimensions with the 
Chern-Simons term becomes tachyonic if we turn on a constant electric field. In contrast, a 
background magnetic field does not cause instability, but it makes the gauge field massive as 
in three dimensions^ 

Chern-Simons terms abound in supergravity theories, and charged black hole solutions 
in these theories provide an interesting laboratory in which to study the instability and its 
implications since these solutions carry background electric fields. The near-horizon geometry 
of the five- dimensional extremal Reissner-Nordstrom black hole in AdS^ is AdSi x M 3 with the 
gauge field strength proportional to the volume form of AdS2- The background electric field 
causes mixing of the gauge field with the metric at the quadratic order, and we will take it into 
account in our stability analysis. We find a critical value a cr ; t of the Chern-Simons coupling 
a above which the near-horizon geometry becomes unstable for some range of momenta k 
in M 3 . Interestingly, the range excludes k — 0, i.e. the instability happens only at non-zero 
spatial momentum. 

1 To our knowledge, [5] is the first paper to point out that the Chern-Simons term in five dimensions induces 
instability. They considered a system consisting of two non-Abelian gauge fields coupled to an adjoint scalar 
field with a tachyonic mass as a holographic model of QCD and reduced it to four dimensions before studying 
its spectrum. Though their set-up and analysis are different from ours, the dispersion relation we derive in 
section [2| is related to theirs. We will point this out at an appropriate place in section [2| 
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Figure 1: Critical temperature as a function of the Chern-Simons coupling a. The shaded 
region indicates a phase with a non-zero expectation value of the conserved current J which 
is helical and position dependent. 

The Reissner-Nordstrom black hole solution in AdS§ gives a holographic description of a 
thermodynamic state in the dual conformal field theory at finite temperature T and chemical 
potential /j^ We find that, for a > a CI a, there is a critical temperature T c (a) below which 
the black hole solution becomes unstable, as shown in Fig. [T] The instability happens at a 
range of momenta, which becomes wider as T is lowered but never includes k = 0, as shown 
in Fig. [2] We find an interesting subtlety in the zero temperature limit; the unstable range 
is wider than the range expected from the analysis near the horizon of the extremal black 
hole. It turns out that the near-horizon analysis gives a sufficient but not necessary condition 
since there are unstable modes in the full Reissner-Nordstrom solution which do not reduce 
to normalizable modes in AdS2 x M 3 in the near-horizon limit. 

In the dual field theory in (3 + 1) dimensions, the instability of the Reissner-Nordstrom 
solution can be interpreted as a signal of a novel phase transition at finite chemical poten- 
tial where the charge current J(x) dual to the gauge field develops a position dependent 
expectation value of the form, 

(J(x)} = Re (ue ikx ) , (1.1) 
with non-zero momentum k. The constant vector u is circularly polarized as, 

k x u = ±i\k\u, (1.2) 

2 El studied the thermodynamic properties of the Reissner-Nordstrom AdS black hole. Its relation to 
Fermi liquid is discussed in [TJ. 
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Figure 2: The left figure indicates unstable regions for various values of the Chern-Simons 
coupling a. The right figure is for a particular choice of the Chern-Simons coupling a = 
1.6a cr it- The critical temperature Tc is the maximum temperature with unstable modes. The 
figure indicates the unstable range a for some temperature T < Tc- The range b is derived 
from the near-horizon analysis at T = 0. Note that the actual range of unstable momenta is 
wider. 



where the sign is correlated to the sign of the Chern-Simons couping as we will explain 



later. The vacuum expectation value (1.1) is helical and breaks translational and rotational 



symmetries in three spatial dimensions, while preserving a certain combination of the two. 
The configuration reminds us of the cholesteric phase of liquid crystals. 

Spatially modulated phases are known in condensed matter physics and in QCD. In the 
Fulde-Ferrell-Larkin-Ovchinnikov phase, a Cooper pair of two species of fermions with dif- 
ferent Fermi momenta condenses with non- vanishing total momentum J9HTU]. An analogous 
effect in QCD was studied in [llj. It has also been shown that finite density QCD in the 
large N c limit is unstable against forming the chiral density wave [121 H3] • Gravity theories 
with Chern-Simons terms may provide dual descriptions of such systems and clarify aspects 
of their phase transitions. 

For example, the Brazovskii model [2] generates a spatially modulated phase, and it has 
been applied to a variety of physical problems (15] . In this model, a non-standard dispersion 
relation is postulated so that the fluctuation spectrum has a minimum at non-zero momentum. 
The gravity theory discussed in this paper provides a holographic realization of a similar 
dispersion relation. 
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In this paper, we use the Maxwell theory with the Chern-Simons term coupled to the 
gravity in AdS§ as a phenomenological model of quantum critical phenomena in the spirit 
of [U |2]. To have an explicit description of the field content and interactions of the dual 
field theory, we need to identify a specific superstring construction where the instability takes 
place. We examined the simplest case of the three-charge black hole in the type IIB superstring 
theory on AdS$ x S 5 and found that the Chern-Simons coupling of the low energy gravity 
theory barely satisfies the stability bound. More specifically, when the three charges are the 
same, the effective Chern-Simons coupling a is only 0.4 % less than the critical value a cr it for 
the instability. There is a limit of an extreme ratio of charges, where an effective a coincides 
with a cr it and the black hole becomes marginally stable. 

This seems to indicate that, if we survey a wider class of examples, we may be able to 
find a theory with a Chern-Simons coupling large enough to cause an instability. Generally 
speaking, the Chern-Simons coupling for a gauge field in AdS$ is proportional to the chiral 
anomaly of the corresponding current in the dual conformal field theory [16]. In particular, for 
the type IIB superstring theory on AdS§ times a toric Sasaki-Einstein manifold, the Chern- 
Simons coupling is determined by the toric data, or equivalently by the combinatorial data of 
the quiver diagram for the dual gauge theory [17J. It would be interesting to find an explicit 
example where the Chern-Simons coupling exceeds the stability bound. Or, one may try to 
prove that such theories are all in the Swampland [HI [19] . 

Even for a theory with a < a cr i t , in which the Chern-Simons term is not strong enough to 
cause instability, the non-standard dispersion relation is noteworthy by itself with potential 
applications to physical problems. For example, a plasmino in QCD is a collective mode of 
quarks whose spectrum has a minimum at a non-zero momentum J20J [21] • In the presence 
of a plasmino, the dilepton production rate diverges at the minimum of the spectrum due to 
the Van Hove singularity, i.e., the divergence of the density of states per unit energy [22] [23]. 

We should also point out that another type of instability of rotating charged black holes 
was suggested in [24"[ 125] . While the Chern-Simons term seems to play a role there, we 
have found no obvious connection to the instability discussed in this paper. Effects of the 
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bulk Chern-Simons terms on hydrodynamics of the dual field theories have been studied in 
[26l [271 1281 [29] . In [301 El]> dispersion relations of hydrodynamic waves in the Reissner- 
Nordstrom geometry with the Chern-Simons term are discussed. Since the authors of these 
papers relied on power series expansions around k — 0, they did not observe the instability 
we found in this paper since the range of instability is away from k = as shown in Fig. 2. 

This paper is organized as follows. In section 2, we show that the five-dimensional Maxwell 
theory with the Chern-Simons term is unstable in the presence of a constant electric field. 
The metric is treated as non-dynamical in this analysis. In section 3, we turn on the metric 
fluctuation and study the stability of the near-horizon geometry of the extremal Reissner- 
Nordstrom black hole in AdS§. In section 4, we generalize the analysis of section 3 to the 
full Reissner-Nordstrom solution. We solve the linearized equations around the black hole 
geometry and identify the critical temperature T cr it of the phase transition. We examine the 
onset of the phase transition and interpret the result from the point of view of the dual field 
theory. In section 5, we show that the three-charge black hole in the type IIB superstring 
theory on AdS 5 x S 5 is barely stable against the type of instability discussed in this paper. 

2 Maxwell Theory with Chern-Simons Term 

It is well-known that the three-dimensional Maxwell theory with the Chern-Simons term is 
massive P, Hj. The equation of motion for the 2-form field strength F is given by 

d*F + aF = 0, (2.1) 

where a is the Chern-Simons coupling constant. Applying d* to this equation and using the 
Bianchi identity dF = 0, one finds 

DF = d*d*F = -ad*F = a 2 F. 

Thus, the Chern-Simons term in three dimensions induces the mass |a| of the gauge field. 

Surprisingly, we find that the Chern-Simons term in five dimensions can turn the Maxwell 
theory tachyonic. In this section, we will demonstrate this by treating gravity as non- 
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dynamical. Coupling to gravity will be studied in the following sections. Consider the 
following Lagrangian density, 

C = -\v^FuF IJ + ^e IJKLM AjF JK F LM , (2.2) 

with the equation of motion, 

-iJI\ i a JJKLM 



dj(^gF J1 ) + —€ FjkFlm = 0. (2.3) 

We use the almost positive convention for the metric gu (I, J — 0, 4). Choose a background 
solution F^ and linearize (2.3) around it by substituting F = F^ + f in (2.3). The linearized 



equation for / is given by 

dj(V=~9f JI ) + ae IJKLM Ff^f LM = 0. (2.4) 



If F^ is magnetic, this equation is similar to (2.1); the fluctuation fjj is massive and the 
configuration is stable. If F^ is electric, on the other hand, (2.4) has tachyonic modes as we 
now explain. 

Suppose the five-dimensional space is flat IR 1 ' 4 , regard it as the product IR 1 ' 1 x IR 3 , and use 
coordinates (x^ =0 ' ,y ' 3 ' ). Let us turn on a constant electric field in the x 1 direction, 

F (0) = Ee 

F^=0, i#> = 0. (2.5) 



The equation of motion (2.4) is then 



d»Uv + d l f iv = 0, 

<)'■])„ + &fji ~ 2aEe ijk f jk = 0. (2.6) 

Our e-symbol convention is such that eoi = 1 and £234 = 1. By multiplying e^dj to the 
second equation, we obtain 

(d^d, + Vdj) fi - AaEt ljk d 3 f k = 0, (2.7) 
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where 



fi ^ijkfjk- 



To derive (2.7), we used the Bianchi identities, 

dif/ij djff^i 9fifiji ^ijkdifjk 29 fi 0. 

In the momentum basis e ip,lX + the operator e^dj has eigenvalues ±k and 0, where 
k = \k\. However, the eigenvalue corresponds to fi ~ k iy which is prohibited by the Bianchi 



identity k % fi = 0. Thus, the linearized equation (2.7) gives the dispersion relation 



(Po) 2 ~ {Pif = k 2 ± AaEk 

= (k ± 2aE) 2 — Aa 2 E 2 . (2.8) 

We find tachyonic modes in R 1 ' 1 in the range of < k < 4 1 aE \ . 

It is instructive to compare this with the case when we turn on a constant magnetic field, 



-(o) _ _ P (o) 

34 
<0) 



-^34 - --^43 - a i 



Ffj = (otherwise). (2.9) 

By repeating the previous analysis, we find the dispersion relation, 

(Po) 2 ~ (Pi? ~ (k 2 ) 2 = (V(h) 2 + (h) 2 + 4a*B* + 2\aB\f . 

In particular, when k 3 = /c 4 = 0, the equation gives pi — p\ — k\ = (4aB) 2 , reproducing the 
topologically massive gauge field in three dimensions. 

In the following sections, we will examine stability of the extremal Reissner-Nordstrom 
black hole in AdS§. If the boundary theory is on K 1 ' 3 , the near-horizon geometry of an 
extremal black hole takes the form AdS^ x R 3 with an electric field proportional to the 
volume form of AdSz- In such a configuration, the effective mass squared in AdSi is again 



3 At this point, we should note that there is a similarity of this dispersion relation to eq. (17) of [5] if we 
set m p = m ai in the paper and interpret m 2 p as being equal to (pi) 2 - 
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given by the right-hand side of (2.8). The configuration is unstable if — 4a E violates the 



Breitenlohner-Freedman bound m BF in AdS2, namely, 



1 



Aa 2 E 2 > \m 2 BF \ = (2.10) 



Ar 2 2 



where r 2 is the curvature radius of AdS^- If this inequality is satisfied, the instability happens 
for non-zero momenta in the range, 



2 i QB i 1 1 - 'Z 1 - i6^kf ) <k< 2 ^ [ 1 + f- } ■ (2 - u » 

It is interesting to note that the zero momentum k = is excluded from the instability range. 
Thus, the condensate of the gauge field happens for non-zero momentum in the M 3 direction 
of the near-horizon geometry. 

As we shall see in the next section, the value of a for the minimal gauged supergravity 



is such that 4a 2 E 2 exceeds the stability bound as in (2.10). This, however, does not mean 
that extremal charged black holes in the minimal gauged supergravity are unstable since we 
must take into account the coupling of the Maxwell field to other degrees of freedom in the 
supergravity theory. We will perform this analysis in the next section. 

3 Coupling to Gravity 

The background electric field causes mixing of the gauge field with the metric at the quadratic 
order, and it modifies the stability condition. In this section, we will study stability of the 
near- horizon geometry of the extremal Reissner-Nordstrom solution in AdS$. It is a solution 
to the Maxwell theory with the Cher n- Simons term coupled to the Einstein gravity with 
negative cosmological constant, 

!QnG 5 £ = V^(r+ 1 4~ ]t 2 FuF IJ ) + ^£ 3 e IJKLM AjFj K F LM . (3.1) 



P 4 1J J 3\ 

The curvature radius r§ of the AdS§ solution in this theory is equal to t. In the following, 
we will work in the unit of I = 1. This is also the Lagrangian density of the minimal gauged 
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supergravity in five dimensions [32]. In this case, supersymmetry determines the Chern- 
Simons coupling a as 

a = —=. (3.2) 
2^3 



In this and next sections, we will treat (3.1) as a phenomenological Lagrangian with a as its 
parameter. 

3.1 AdS 2 x R 3 

Let us first consider the extremal black hole solution which is asymptotic to AdS$ in the 
Poincare coordinates. It describes the dual conformal field theory on IR 1,3 with non-zero 
chemical potential and at zero temperature. The near-horizon geometry of the extremal 
black hole is AdS2 x IR 3 with the metric 



ds = ,n/ n, + d f, y = (y\y d y). (3.3) 



-jd x ) 2 + (d x 1 ) 2 

nix 1 ) 2 

Note that the curvature radius of AdS2 is 1/ vT2; the curvature is stronger near the horizon. 
The electric field strength near the horizon is proportional to the volume form of AdS% and 
is given by 



For the minimal gauged supergravity with a given by (3.2) 



Aa 2 E 2 = 8 > \m 2 BF \ = = 3. (3.5) 

Thus, if gravity is treated as non-dynamical, the gauge field fluctuation near the horizon 
violates the Breitenlohner-Freedman bound for this value of a. 

We decompose the metric gjj into the background gf] and the fluctuation hu as gjj = 
9lj + hi J- The indices are raised/lowered by using the background metric. Notice that 
g IJ = g IJ (°) — h IJ + 0(h 2 ) so that g IJ gjK — ^k- m ^he presence of the background electric 
field Ff£/ , the unstable gauge field components f^, ^ mix with the off-diagonal elements 
h 2 of the metric perturbation through the gauge kinetic term, 

F U F IJ = 4F^hjf ui + ■■■ . (3.6) 
9 



Thus, in the stability analysis, we have to take into account the mixing. One can think of 
as the Kaluza- Klein gauge field upon reduction on M 3 . Since we are considering a sector with 
non-zero momentum k along IR 3 , the Kaluza-Klein gauge field on AdS 2 has mass k 2 . 

To examine the stability of the black hole solution, we can apply the standard linear 
perturbation theory. In the present situation, however, there is a simpler way as we describe 
here. Suppose that the momentum k on M. 3 is in the y 2 direction. To derive the effective 
action for the Kaluza-Klein gauge field in AdS2, it is convenient to reduce the Einstein 
action in (3.1) along the y 3,4 directions first. This gives rise to two gauge fields {hj,h 2 l ) 
(i = 3,4) on AdS 2 x M. y 2, with the effective Lagrangian 



atl (R + 12) 



(0) 

-9sd 



i=3,4 



(^K^K^" + ^K^K^ + (terms not involving hj, hj) 



(3.7) 



where the gauge field strengths are 



k; v = d,h; - d„h;, k; 2 = - d 2 h; ^ v = o, i ; < = 3, 4). 

Upon further reduction in the y 2 direction with momentum k, the effective Lagrangian density 
for the Kaluza-Klein gauge field is 

"1 



eff 



i=3,4 



(3.8) 



We see that the off-diagonal elements (i = 3, 4) of the metric fluctuation give rise to two 
massive gauge fields of mass \k\ on AdS 2 with h 2 l serving as the requisite Stiickelberg fields. 

Let us dualize the Kaluza-Klein field strength K l on AdS 2 and write it as a function Ki 
times the volume form, 



K, 



01 



I2(x 



The equations of motion for /j = \eijkfjk and Kj are derived from the Lagrangian density 



which is (2.2) plus (3.8) with the coupling (3.6). They can be organized into the form 



(D AdS2 + d J dj) fi - AaEeijkdjfk + Ee ijk djK k = 0, 
E n AdS Ji + (n AdS2 + d 3 dj) e ijk djK k = 0. 



(3.9) 
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The effective mass m of these fields in AdS-z can then be computed by solving 

/ m 2 - k 2 - 4aEk E \ 

det { Em 2 m 2 -k 2 ) = °> ^ 

where k = ±|/c|. We find 

^ (2k 2 + E 2 + AaEk ± v/E 4 + 8aE*k + 4(1 + Aa 2 )E 2 k 2 ^ . (3.11) 



m 2 , 



Minimizing m 2 with respect to k and choosing the minus sign in (3.11), we obtain the 
lowest value of m 2 as 



m 2 



E 2 (-64a 6 - 24a 4 + 6a 2 - (16a 4 + 4a 2 + if' 2 + 1 



2(4a 2 + l) 2 

Substituting E = 2y / 6 for the near-horizon geometry, we find numerically that the lowest 
value of m 2 violates the Breitenlohner-Freedman bound if 

|a| > a crit = 0.2896 ••• . (3.12) 

The value of a for the minimal gauged supergravity is 

a = — = = 0.2887- ■ ■ . 

Thus, the supergravity theory is stable against the fluctuation of the gauge field, but barely 
so (with a margin less than 0.4 %). 

3.2 AdS 2 x S 3 

For completeness, let us consider the case when the boundary theory is on K x S 3 . The 
near- horizon geometry is AdS% x S 3 . Let us denote the curvature radii of AdS2 and S 3 by T2 
and rs, respectively. They are related to the electric field strength E and the cosmological 
constant A, which is —6 in the AdS2 x H^ 3 limit, by 

1 2 

A = ~2? + 7 2 ' 
zr 2 '3 

£ 2 = A + A (3.13) 

r 2 r 2 
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Note that, in the limit r 3 — > oo, where S 3 becomes M 3 , this reproduces E = ±2v^6 in our 
unit. 

As in the previous case, we consider fluctuations of the metric gjj = gf] + hjj and the 

? (o) 
u 

Einstein equation, 



gauge field Fjj — Ff) + fjj from their classical values indicated by We expand the 



Ru - \guR = \ (f ik Fj k - - A g u F KL F KL ) . ( 3. 1 i ) 



and the Maxwell equation modified by the Chern-Simons term, 

v^V jF JI + ^e IJKLM F JK F LM = 0, (3.15) 

to the linear order in hu and fu. 

The linearized equations for and K i: where Ki is defined such that i^j(vol AdS2) Mi / = 
2V[ M /i i/ ]j, can be written as, 

(n AdS2 + A 5 a) / - 4aEd*f + EdK = 0, 

ED AdS J + (n AdS2 + A 53 + -|j dK = 0. (3.16) 



These equations are similar to (3.9), except for the last term -j in the second equation. Here 

' r 3 

* means the Hodge dual on S 3 . Since d* is hermitian when acting on the space of two-forms 
on S 3 , decompose / into its eigenstate. Its eigenvalue is known to be k = ±(n + 2)/r^, where 
n = 0, 1,2, .... Since A S 3 = —(d*) 2 when acting on / satisfying the Bianchi identify df = 0, 
we can set A S 3 = —k 2 . 

The mass m on AdS2 then satisfies the determinant equation, 

(m 2 -k 2 -AaEk E \ , . 

d6t ^ Em 2 m 2 -k 2 + A/r 2 J ~ U - [6 ' U) 

This can be solved to obtain, 



2 1 

m = - 



2k 2 + E 2 + AaEk - ^ ± \jE A + 8aE 3 k + ™ + + AE 2 ( k 2 + 4a 2 k 2 - \ 



(3.18) 



In the limit of r 3 — > oo, this reduces to the previous result (3.11) 
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We have numerically checked that, for a wide range of A and E, the Breitenlohner- 
Freedman bound is not violated in the minimal gauged supergravity, where a = It 



is interesting to note that, in the limit of A — > but with non-zero E, the lowest m 2 in (3.18) 
saturates the Breitenlohner-Freedman bound [33], which is 

1 1 / „o . v E 2 

-¥T-I2< E - 2A > = -T2- (3 ' 19 > 

4 Phase Transition and Critical Temperature 

In the last section, we studied the instability of the near-horizon region of the extremal 
Reissner-Nordstrom solution. This gives a sufficient condition for the solution to be unstable. 
However, as we will see in this section, the condition turns out to be not necessary. To 
clarify the nature of the phase transition and identify the critical temperature, we study 
linear perturbation to the full Reissner-Nordstrom black hole in AdS^,. 

4.1 Geometry and Equations 

The Reissner-Nordstrom black hole has the metric 

ds 2 = -H{r)dt 2 + -^-dr 2 + r 2 dy 2 , y = (y 2 ,y 3 ,y 4 ) . (4.1) 
H (r) 

Note y—g^ = r 3 . The gauge field strength is given by 

= Qdt A dr . (4.2) 



The function H(r) is given by 



H(r) = r 2 



1-11 + ^ ( r± Y+ /r+x " 



3ri / \ r / 3ri V r 



(4.3) 



where Q = — 2 / ur 2 h . 

The equation of motion coming from the variation of the gauge field is 

d,(V^r) + ^(v^f) - 2a^e idk f jh - d p (^/^W^h^ = . 



(4.4) 
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In the black hole background (4.1), it becomes 



' d t f ti + d r (rH(r)f ri ) + l-djfji - 2a^e ijk f jk + QK' 1 = 



H r 



where K % = d t h r l — d r h t % . By operating e^dj on this equation, we obtain 

- j^ 9 tfi + d r (H(r)rd r fi) + ^Ans/i - \o ( ',<,,,<),]),. + Qe ijk djK k = 



(4.5) 



(4.6) 



where / { = \e ijk f jk and A M 3 = d 2 2 + d 2 3 + d 2 4 . 

To obtain the equation of motion that comes from the variation of the off-diagonal metric 
elements, let us use the Kaluza-Klein reduction in the presence of momentum k along the y 2 
direction. The effective Lagrangian has the form 



C 



eff 



—r 



i=3,4 



\r 2 K\K iiXV + - \djij - ikh*\ 2 



(4.7) 



The r 3 factor comes from the volume form on the M 3 directions with coordinates y. The 
equation of motion coming from the variation with respect to the metric is given by 



d u (r 5 K u ^) + r*(-ik)(d ft h 2 i ~ ikh ^) ~ Q^ P f P i = • 



Acting on the operator ea^dpg" 13 ^^, we can eliminate the term containing d fM h 2 t . Using 
+ e^da + e va d^ = in two dimensions, we obtain 



- ^dpg^^d^K^ + k 2 K l - Qd^ ^^) = . 
Further operating e^dj on the equation, 

dp[g am ^d a {rh ijk d,K k )} + A 2 R3 e ijk d 3 K k + Qd^g^^dJ,) = . 
More explicitly, 

■^rft + d r H(r)±d r )j {r% jk d K k + Qf t ) + A^e^tf* = . 



(4.9) 



(4.10) 



(4.11) 



We have two sets of equations of motion (4.6) and (4.11). To simplify them, let us perform 
the following rescaling, 



r+ 
u 



14 



y 



X 



(4.12) 



and make the change of variables, 

fi(r) -> 0(r), 



e ijk djK k -> ^-«V(r), (4.13) 



and set g = . The temperature T is 



r=£(2-£)- (lid 



2tt V 3r 

With the rescaled variables, the Reissner-Nordstrom black hole is 



.2 ' ( //, ..»./, 2 , 1 J„.2 \ , 1 j -2 



ds 2 = — -H(u)dt 2 + -= du 2 + —dx 2 , (4.15) 

u V H(u) J u 

where 

ff(ti) = i- (i + g y+?v. (4.i6) 

In these coordinates, the AdS$ boundary is at u = and the black hole horizon is located at 
u = 1. 

Suppose that the fields and ^ have time dependence e~ wt . Then the equations of motion 
for the fields and ip give the following set of ordinary differential equations, 



cu 2 



H(u) 



+ ud u [H{u)u~ l d u (^j - k 2 (j) + 8V3aqku 2 (f) - 2qu 2 ^ = 



CO 



2 



(ip - Qqu 2 (f)) + u~ 1 d u H(u)u 3 d u (u- 2 (iP - Qqu 2 <f))) 



(4.17) 



-k 2 ij = o 



H(u r 

Introducing a new function £ = ip — 6qu 2 <f>, the equations can be written as 
+ ud u (H(u)u- l d u (p) - (k 2 + 8VSaqku 2 + 12gV)0 - 2qu 2 £ = 



cu 2 



Hun 

-f + w9 u (#(w)w _1 9 u - 6qk 2 u 2 <j> - (A; 2 - 8m 2 - 9u 2 q 2 + 12gV)£ = . 



a; 2 



(4.18) 



H{u) 

Interestingly, the two equations can be diagonalized by a -u-independent matrix. That is, for 
some linear combinations 0i an d 02 of and £, we have 

w 2 

7 r— 0i(w) + ud u (H{u)u~ l d u (t)i{u)) - K,i(u)(j)i(u) = , (4.19) 
if(u) 
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where % — 1 , 2 and 

Ki(u) =k 2 — 4\/3akqu 2 

- 2u 2 ^2 + g 2 (2 - 6m 2 ) =f \A + 4g 4 - 8y/3akq - 8V3kq 3 a + q 2 (S + P(3 + 12a 2 ))^ , 

(4.20) 

where k\(k2) chooses the minus (plus) sign on the right-hand side. Our numerical analysis 
shows that only 02(w) can be an unstable mode. It is related to the fact that, in the extremal 



limit, K2 corresponds to the smaller mass-squared in (3.11) in the near-horizon limit 



4.2 Numerical analysis 



To solve the equations of motion (4.19) numerically, we impose the in-going boundary condi- 



tion near the horizon u — 1, and then evolve the solution to u — 0, the AdS§ boundary. The 
asymptotic behavior of <pi near u = is either 0j ~ u 2 or constant. The former is normalizable 
and the latter is non-normalizable. To find normalizable modes in the full Reissner-Nordstrom 
solution, we scan the initial conditions and see when the fields vanish at u — 0. 




■ a/a clit 



Figure 3: For a given value of the Chern-Simons coupling a, there is a discrete set of momenta 
k for which static solutions exist. The curves I and II indicate two of such momenta for each 
a. The red curve is the lower-end of the momentum range that violates the Breitenlohner- 
Freedman bound near the horizon. Note that the red curve coincides with the curve II. 
However, there is another curve I with a lower momentum. This means that the near-horizon 
analysis gives a sufficient but not necessary condition for the instability. Both curves end at 
the same critical value of a. 



First, let us consider the zero temperature limit (q = \/2) and search for static solutions 
{u = 0), which signal the onset of an instability. The behavior of the fields near u = 1 can 
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be found from (4.19) as 

0i = (l-«)-3+V ^^(1 + •••), (4.21) 

where terms in • • • vanish at u — 1. In the actual numerical calculation in this section, we 
include several subleading terms to improve accuracy. For a given Chern-Simons coupling a, 
static modes appear at discrete values of momentum k. The lowest two modes are plotted in 
Fig. |3j As mentioned before, only the second field 02 has normalizable static solutions. 

The two curves in Fig. [3] are denoted as I and II. Both curves terminate at a/a CV i t = 1 
and k/fi = 1.52 • • • . The critical value of the Chern-Simons coupling a cr i t = 0.2896 • • • is the 
one we found from the stability analysis of the near-horizon geometry in the previous section. 
The curves are supposed to extend over k/(i = 1.52- •■ and come back to the right in a 
bell-shaped curves. The upper branches of the curves represent the upper bounds of unstable 
modes. However, we have not been able to plot them due to inaccuracy of our numerical 
computation. 

We also found out a static solution at zero momentum. However, for this solution, the 
curl of the off-diagonal metric component €ijkdjK k is constant on M. 3 . This means that K % is 
linear in M 3 , and the solution is not normalizable. 

We note that the curve II fits with the red curve which is at the lower-end of the momen- 
tum range that violates the Breitenlohner-Freedman bound in the near-horizon AdSi x M 3 
geometry. As we saw in the previous section, the near-horizon geometry is unstable in this mo- 
mentum range, thus the full Reissner-Nordstrom solution should also be unstable. In fact the 
momentum range that violates the Breitenlohner-Freedman bound is specified by k 2 (1) < —3, 
where §i{u) oscillates infinitely many times as they approach the horizon as can be seen from 



(4.21). On general ground, we expect an instability to occur in this range 

Interestingly, the instability condition /^(l) < — 3 of the near- horizon geometry is not 
necessary for the instability of the full solution. This is because there is yet another curve I, 
located outside of this momentum range. What happens is that the curve I corresponds to 
a normalizable perturbation to the full Reissner-Nordstrom geometry, but the corresponding 
mode becomes non-normalizable in the near-horizon limit. Our numerical analysis shows that 

17 



the critical Chern-Simons coupling a cr i t for the curve I is the same as that for the curve II, 
even though the value of a cr it was derived from the near-horizon analysis. 

To see that these static solutions indeed signal instability, let us turn on u with positive 



imaginary part in (4.19). We impose the in-going boundary condition, which is 



e 12 < 



^)(l-u)^ M (i + 



(4.22) 



in the zero temperature limit and 



(4.23) 



at a positive temperature. Fig. [| shows the negative frequency squared as a function of the 
momentum k at zero and finite temperature. It shows that the upper and lower curves in 
Fig. [3] are boundaries of unstable modes. 

T=Q -10 8 wV r = 8.7xHr'/' 



- 10V//i 
k 

3 - 

2 - 



- 10V//i 
A 



k/fi 



k/fi 



Figure 4: Left: Negative frequency squared as a function of momentum k at zero temperature 
when a = 1.6a cr ; t . Only positive — u 2 is plotted. The curves starting around 1 and 3 join to 
represent a tachyonic dispersion relation for the unstable mode predicted by the near-horizon 
analysis. The curve starting below 1 is also expected to be connected with another line in 
the higher momentum region to form a larger bell-shaped curve, but the large momentum 
part is difficult to analyze numerically. The zero momentum static solution does not extend 
to an unstable mode. Right: Negative frequency squared as a function of momentum at 
temperature T = 8.7 x 10~ 4 /z- 



The occurrence of instability by the Chern-Simons coupling is summarized concisely in 
Fig. [T] and Fig. [2] in the introduction section of this paper. For each Chern-Simon coupling 
a, Fig. [2] shows an unstable region in the momentum-temperature plane. This is related to 
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the curve I in Fig. [3] The range of unstable momenta never includes k = 0. The highest 
temperature with unstable modes is denoted as Tc(a). Fig. [I] shows this critical temperature 
as a function of a. Below the critical temperature Tc(oc), we expect an instability and the 



charge current gets a position dependent expectation value of the form of (1.1). 
4.3 Spontaneous Current Generation 

The vacuum expectation value of the current J in the dual field theory can be evaluated 
by extracting the asymptotic behavior of the corresponding gauge field toward the boundary 
of AdS$. In the absence of the Chern-Simons term, it is well-known that (J) is given by 
the normalizable part of y—g^f n evaluated at r — ► oo. The normalizable mode of / r , 
decays, but that effect is compensated by the scaling behaving of the metric so that we find 
a finite limiting value in the low temperature phase. The Chern-Simons term gives rise to 
an additional term of the form afieijkfjk- However, it vanishes at the boundary and does not 
contribute to the expectation value. Thus, the vacuum expectation value of the current in the 



low temperature phase is given by ^/—g^f Tl evaluated at the boundary of AdS§. It takes 
the form, 

(J(x)) = Re (ue ikx ) , (4.24) 
where the polarization vector u obeys 

k x u = ±i\k\u. (4.25) 

From the analysis in the previous section, it is clear that we should choose the plus (minus) sign 
when aE is positive (negative). Namely, the sign of the Chern-Simons coupling determines 
whether the circular polarization of the current expectation is clockwise or counter-clockwise. 
This configuration breaks translational and rotational symmetries, but a certain combination 
of the two is preserved. The polarization of the current is helical and reminds us of the 
cholesteric phase of liquid crystals. 

Since the gauge field mixes with the metric fluctuation in the bulk, the correspond- 
ing component T 0i of the energy-momentum tensor has a non-zero expectation value at the 
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boundary. This is expected since the non-zero current in the spatial direction means that 
there is a momentum density. 

4.4 Spontaneous Breaking of Internal Symmetry 

So far, we have considered the case when the gauge group in the bulk is U(l). Since the U(l) 
current commutes with itself, its expectation value does not break the U(l) global symmetry 
on the boundary. 

To realize spontaneous breaking of an internal symmetry, one possibility would be to 
choose the gauge group to be non-abelian. The Chern-Simons term can be written in five 
dimensions if there is a symmetric tensor d a b c in the Lie algebra, such as in 577 (n) with n > 3. 
Suppose we turn on an electric field strength in a direction T a in the Lie algebra. According 
to EE] , the gauge kinetic term can generate instability in directions in the Lie algebra that 
do not commute with T a . This breaks the symmetry homogeneously. On the other hand, the 
Chern-Simons term can cause a spatially modulated instability in directions where d a j, c ^ 
with T a . The competition of the two effects would be decided by the relative strength of 
the gauge coupling and the Chern-Simons coupling. It would be interesting to study such an 
effect in a more explicit manner to identify the gravity dual of a spatially modulated phase 
with spontaneous breaking of an internal symmetry. 

5 Three- Charge Black Holes in Type IIB Theory 

The consistent truncation of the type IIB theory on AdS 5 x S 5 to the U(l) 3 gauged super- 
gravity in five dimensions was given in [37J . The bosonic action contains three gauge fields for 
£7(1) 3 and three scalar fields Xi, X 2 , X 3 subject to the constraint XiX 2 X 3 = 1, in addition to 
the metric. This low energy theory admits the three-charge black hole solutions of [38]. Here 
we will examine the stability of the near-horizon region of the three-charge black holes in the 
extremal limit. 
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5.1 Case with Equal Charges 

Let us consider the case when the three charges are identical, which implies the scalar fields 
are constant X 1 — X 2 — X 3 — 1. In this case, both the Lagrangian and the black hole con- 
figuration are symmetric under exchange of the three gauge fields Fi,F 2 ,F 3 . It is convenient 
to take their linear combinations as 

F = ~^=(F 1 + F 2 + F 3 ) 

F+ = ^= {F 1 + uoF 2 + u 2 F 3 ) 
v3 

F_ = 4= (*i + u 2 F 2 + uF 3 ) . (5.1) 

They are eigenstates of the Z 3 permutation with eigenvalues 1 and u^ 1 , where u = e 2nl ^ 3 . 
In the black hole geometry, the Z 3 invariant gauge field F has an electric component with 
E = 2\/6, and F± — 0. Similarly, fluctuations of the scalar fields from Xi = X 2 = X 3 = 1 
can be organized into eigenstates with eigenvalues uj ±x under the Z 3 permutation. 

The Z 3 invariant sector is the minimal gauged supergravity with a = l/2\/3. To the 
quadratic order, the metric and the Z 3 invariant gauge field do not mix with other fields. 
Thus, the stability analysis with respect to them is exactly the same as the one we performed 
in the previous section. The three-charge black hole is barely stable in this sector, being 
within 0.4 % of the stability bound. 

Since the gauge fields F± have zero expectation value on the black hole geometry, the M 3 
components of these gauge fields do not couple with other degrees of freedom in the quadratic 
order. It is convenient to write them as 

F ± = i=(/«±^ 2 )). 

With the standard normalization of their kinetic terms, the Chern-Simons term takes the 
form, 

1 TJKLMjp / (l) f (l) , (2) f (2)\ /- 9 x 

e fiu a K Jlm + a K Jlm I > \P- Z ) 



8^/3 

where are the vector potentials for ffj (i = 1.2). To the quadratic order, we can take the 
Z 3 invariant Fjj to be its background value Fjj . 
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Since these gauge fields do not couple to other fields in the quadratic order, their linearized 
equations of motion are, 

dj(^9f (i)JI ) + ^ JKML Ffi& = 0, (a = 1, 2). (5.3) 



Comparing this with ( 2.4 ) , we find a = ±1 /Ay/3. Since as in the previous example, 

the mass squared is given by 

- Aa 2 E 2 = -2. (5.4) 

It is greater than the Breitenlohner-Freedman bound (m 2 BF = —3 in our unit), and quadratic 
fluctuations in these gauge fields are stable. 

5.2 Case with Non-Equal Charges 

Next, let us consider the case when the three charges are different. The five-dimensional 
Lagrangian is derived in |37j . 



16ttG r £ 



\ a a / 



JJKLMttiI 7ti2 A 3 
4 e ^IJ^KL A M ■ 



(5.5) 

X a are functions of the two scalars (f)% and 4> 2 subject to the constraint XiX 2 X 3 = 1. The 
Lagrangian admits AdS$ black holes parametrized by three charges qi, q 2 and q%. The metric 
is given by 

dr 2 



i 



ds 2 = -{H x H 2 Hz)-ih{r)di 2 + (HiH 2 H 3 )^ ( ^ + r W 

la . ,2 



Mr) ' ' 3 
H a (r) = l + ^, q a = H sinh z (3 a , a = 1,2, 3 

X a = H-\H l H 2 Hz) 1 * 

h{r) = l-^- + r 2 H x H 2 H^ 



r 



A, = (l-^'r 1 )coth^ a tft. (5.6) 
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This is the metric whose foliating transverse space is S 3 . If it is IR 3 instead, the S 3 metric 
dfll is replaced by the flat metric and h(r) and A a are replaced with 



h(r) = -^+r 2 H 1 H 2 H 3 



I -HZ 



-dt . 



(5.7) 



sinh/3 a 

Given the charges q a , it may be possible to choose p such that the black hole becomes extremal. 
That is, the inner and the outer horizons coincide. For the extremal case, the near-horizon 
geometry is AdS 2 x S* 3 or AdS 2 x IR 3 : if the horizon occurs at r = r , for small p = r — r , 
h(r) = lh"(r )p 2 . Hence the geometry becomes 



ds 2 



1 

<2l 



2 >,2 . 

-p dt H 7T- 

p 2 



H dQ 2 3 , 

a 2 



(5.8) 



where a\ = ^(HiH 2 H 3 )~^h" , and a2 = (HiH 2 H 3 )~3rZ~ 2 for S* 3 and is replaced with 

the flat metric for IR 3 . H a and /i" are implicitly evaluated at r = r . 

We want to analyze the linear fluctuations near the horizon in the extremal limit. Let 
F a = Fa°^ + f a = Fa°^ + da a and gu = + hjj. If we focus on the fluctuations of the Oj and 

fields only, we find that the linear fluctuations of the scalar fields 0i and (f> 2 do not couple 
to them. Therefore, we may use the background value of the scalar fields. In this case, we 
can derive the equations of motion as in the previous case, and the result is 

(A 2 + A3)/ 1 - X 2 E 3 d*f 2 - X 2 E 2 d*f 3 + E 1 dK = 
- X 2 E 3 d*f + (A 2 + A 3 )/ 2 - X 2 E ± d*f 3 + E 2 dK = 



- X 2 E 2 d*f L - XiE 1 d*f 2 + (A 2 + A3)/ 3 + E 3 dK = 
|^A 2 7 1 + ||A 2 7 2 + ||A 2 7 3 + (A 2 + A 3 + Aa 2 )dK = . 

A l A 2 A 3 



(5.9) 



E a are the electric fields such that dA a = E a (HiH 2 H 3 ) sdt A dr. The above four equations 
give a mass matrix equation 



det 



(m 2 -k 2 -E 3 X 2 k -E 2 X 2 k 
-E 3 X 2 k m 2 -k 2 -E x X 2 k 



Ex 
E 2 
E 3 



\ 



-E 2 X 2 k -E x X 2 k m 2 -k 2 
\ fj m2 % m2 f| m2 m 2 -k 2 + 4a 2 J 



= . 



(5.10) 
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Solving this equation for m 2 , we obtain the mass spectrum. 

When two of the three charges are the same, we can analyze the mass spectrum analytically 
for the AdS 2 x M 3 geometry. In this case, only the ratio of the charges matter. Let the charge 
assignments be (q a ) = (l,q,q). Demanding /(r ) = f'(r ) = at some r = r , we obtain the 
relation q = x(2x + 1) and \i = Ax (1 + a;) 3 where x = r\. Let us parametrize the extremal 
solutions in terms of x. Then the various functions at the horizon are given by 

H x = X -\l + x) , H 2 = H 3 = 2(1 + x) , H 1 H 2 H 3 = Ax~\l + x) 3 (5.11) 

X x = 2M , X 2 = X 3 = 2"52T5 (5.12) 

E x = 25x5, E 2 = E 3 = 23x"5(l + 2x)5 (5.13) 

a x = 23x"s(l + Ax) , A = -2sx"3(l + Ax) . (5.14) 

When the two charges are the same, there is a Z 2 symmetry exchanging the two charges. 
Since the gravity is insensitive to this exchange, only the combination f 2 + f 3 couples to the 
metric component and f 2 — f 3 decouples. The decoupled mode is analyzed by considering 



the mass matrix in (5.10) with the eigenvector (0, 1, —1, 0) for some m 2 . Due to the fact that 
E 2 = E 3 and X 2 = X 3 , the only condition that we need to satisfy is 

m 2 - k 2 + E x X\k = . (5.15) 

Therefore m 2 has the minimum value when k = El ^ 2 , in which case m 2 = — El *' 2 = —2^x^. 
The Breitenlohner-Freedman bound is — ^ = — 2 _ ix~i(l + Ax). Their ratio is jt|^, which is 
always lower than 1. That is, the mass squared is always above the bound. 

Of course, it is possible that there are other modes that go below the bound. But this 



turns out not to be the case. To see this, let us evaluate the determinant (5.10) when the 
mass-squared m 2 takes the value — \a\, which is the Breitenlohner-Freedman bound. Then 
this is a function of k and x. We can check that this function is always positive, meaning that 
the roots of the determinant equation, which are the possible values of the mass-squared, are 
all greater than the Breitenlohner-Freedman bound. 

When all three charges are different, we have not been able to solve the equations an- 
alytically, so we resorted to a numerical method. Given three charges, we first adjust the 
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parameter /i in (5.6) so that it gives an extremal black hole. Then we evaluate the metric 
and the functions at the horizon and solve the mass matrix equation (5.10) for m 2 . In both 



AdS% x S 3 or AdS% x M 3 , however, no unstable modes are found for a large range of the three 
charges. The bound is always barely satisfied. 

In this section, we studied stability of the three-charge black hole in the near-horizon 
limit. As we saw in the previous section, the near-horizon analysis gives a sufficient but not 
necessary condition for the instability at T = 0. However, the critical value of the Chern- 
Simons coupling is given correctly from the near-horizon analysis. Thus, we expect that our 
conclusion in this section would not be modified even if we perform the analysis in the full 
black hole geometry. 
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